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The influence of molecular hydrogen on propene polymerization promoted by the regio- 
and isospecific rac-Me 2 C(3-f-Bu-l-Ind)2ZrCl 2 /MAO catalytic system has been investigated 
both by experiment (liquid propene, T p = 50 °C) and molecular modeling (quantum 
mechanics/molecular mechanics). A competitive secondary (2,1) propene insertion into the 
Zr-H bond, originated after chain transfer to hydrogen, is detected by 13 C NMR analysis of 
these polymers: increasing amounts of the new CH3— CH(CH3)— CH(CHa)— CH2— chain-end 
group are observed when the hydrogen partial pressure is increased from 1 to 4 bar. The 
ratio between primary (1,2) and secondary initiations is found to be constant and nearly 
equal to 5. On the other hand, the complete regiospecificity of propene insertion during chain 
propagation is preserved. A combined quantum mechanics/molecular mechanics study on 
the model cation [rac-H 2 C(3-f-Bu-l-Ind) 2 Zr-H(propene)] + shows that secondary propene 
coordination, followed by secondary insertion into the Zr-H bond which generates the 
observed structure of the saturated end group, is even favored over primary coordination/ 
insertion, thanks to a stabilizing /?-agostic bond between the metal center and one methyl 
hydrogen of the secondary coordinated propene. Two different mechanisms have been 
identified for primary propene insertion into the Zr(/>propyl) and Zr(isopropyl) species, with 
energy barriers of nearly 4 and 9 kcal/mol, respectively. The experimentally observed ratio 
of [1,2]/[2,1] ^ 5 for initiations is explained in terms of side reactions that can transform 
the "slow" propagating Zr(isopropyl) site into the "fast" Zr(/>propyl) one. The model [rac- 
H 2 C(3-f-Bu-l-Ind) 2 Zr(isopropyl)] + cation shows indeed comparable energy barriers for 
monomer insertion and £-H transfer followed by associative displacement with a primary 
propene to generate the "fast" propagating species. 



Introduction 

Molecular hydrogen is the most used chain transfer 
agent for the control of polymer molecular weights in 
Ziegler-Natta polymerization of 1 -olefins. 23 Chain trans- 
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(2) (a) Guastalla, G.; Giannini, U. Makromol Chem., Rapid Com- 
mun. 1983, 4, 519. (b) Hayashi, T.; Inoue, Y.; Chujo, R.; Asakura. T. 
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fer most likely occurs via direct addition of H2 across 
the metal-polymer bond with formation of a saturated 
polymer chain and a metal— hydride bond: 4a ~ e 

L n MR + + H 2 — ► — LnMH* +RH 

L n M~ R 

Besides lowering the molecular weight of the polymer, 
addition of hydrogen generally causes an increase of 
polymerization activity using both heterogeneous 2 (MgCV 
TiCl 4 ) and homogeneous 3 (metallocene) catalysts. The 
most common explanation of this activating effect is the 
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CarwiU, A. Tritto, I. Locatelli, P.. Sacchi, M. C. Macromolecules 1997, 
30, 7056. 

(4) (a) Brintzinger, H.-H. J. Organomet. Chem. 1979, 171, 337. (b) 
Gell, K. I.; Posin, B.; Schwartz, J.; Williams, G. J. Am. Chem. Soc. 
1982, 104, 1846. (c) Wochner, F.; Brintzinger, H.-H. J. Organomet. 
Chem. 1986, 309, 65. (d) Guo, Z.; Bradley. P. K.; Jordan, R. F. 
Organometallics 1992, 11, 2690. (e) Ziegler, T.; Folga. E.; Berces, A. 
J. Am. Chem. Soc. 1993. 115, 636, and references therein. 
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3 Polymerization conditions: 1 L stainless steel autoclave, 
propene (0.4 L), 50 °C, 1 h, 0.5 mgof rac- [Me 2 C(3^Bu-Ind) 2 ]ZrCl 2 , 
solution zirconocene/MAO aged 10 min, Al/Zr = 5000 mol/mol. 
^Hydrogen partial pressure, In THN, 135 °C. ^From [rj\ = K(M v ) a 
with K= 1.93 x 1(T 4 and a = 0.74. 'From 13 C NMR (100 MHz). 
^Determined assuming the enantiomorphic site model. 

regeneration of an active species from the "dormant 
sites" formed after a regioirregular olefin insertion. 2 3 - 5 
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slow propagating 
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fast propagating 
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In fact, 13 C NMR spectra of polypropene samples 
prepared in the presence of hydrogen with Ziegler- 
Natta catalysts reveal the presence of a greater amount 
of n-butyl groups 2 ^ 3 *' 5 (generated from chain transfer 
after a secondary insertion) with respect to isobutyl 
groups (from chain transfer after a primary insertion) 
as terminating polymer chain-end. The initiating poly- 
mer chain-end invariably present in these spectra is 
n-propyl which is due to the 1,2 propene insertion in 
the metal-H bond. 

We recently reported on the polymerization of propene 
with the rac-Me 2 C(3-rerf-butyl-l-indenyl)2ZrCl 2 /MAO 
(1/MAO) catalytic system, 6 which produces fully regio- 
regular, highly isotactic polypropene (e.g., mmmm = 
95%, T m = 152 °C at T p = 50 °C in liquid propene). 

In this work we discuss the influence of hydrogen 
addition on liquid propene polymerization with 1/MAO. 
The experimental results are compared with those from 
a combined quantum mechanics/molecular mechanics 
study on insertion and isomerization reactions. 

Results and Discussion 

Propene Polymerization and Analysis of Satu- 
rated End Groups. In Table 1 we report the polym- 
erization conditions and the relevant results obtained 
in the polymerization of propene with rac-Me 2 C(3-f-Bu- 

(5) Busico. V.; Cipullo, R.; Corradini, P. Makromol Chem. Rapid 
Commun. 1992, 13, 15. 

(6) (a) Resconi, L.; Piemontesi, F.; Camurati, I.; Sudmeijer, O.; 
Nifant'ev, I. E.; Ivchenko, P. V.; Kuz'mina, L. G. J. Am. Chem. Soc. 
1998. 120, 2308. (b) Resconi, L.; Polym. Mater. Set Eng. 1999, 80, 421. 

(c) Resconi, L.; Camurati, I.; Sudmeijer, O. Top. Catal. 1999. 7, 145. 

(d) Camurati, I.; Fait, A.; Piemontesi. F.; Resconi. L.; Tartarini, S. In 
Advanced Catalysis: New Polymer Syntheses and Modifications, Boffa, 
L. S., Novak, B. M., Eds.; Symposium Series; American Chemical 
Society: Washington, DC, in press. 
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Figure 1. (a) 13 C NMR spectra of polypropenes prepared 
with Me 2 C(3-f-Bu-Ind) 2 ZrCl 2 /MAO: I = sample 1, II = 
sample 6. End groups: A = r?-propyl, B = isobutyl, C = 
2,3-dimethylbutyl. (b) Expanded methyl pentad region of 
sample 6. 

Ind)2ZrCl2/MAO in the presence of different amounts 
of hydrogen. 

As normally observed, the polymer molecular weight 
is lowered by adding hydrogen. Furthermore, the ad- 
dition of hydrogen produces up to a 2-fold increase in 
productivity. 6d This differs from literature data on 
heterogeneous 2 systems, where a stronger activating 
effect is generally reported, but is in line with what is 
observed for homogeneous, less regiospecific catalytic 
systems. 3 The fact that H 2 activates both highly re- 
giospecific and less regiospecific zirconocene catalysts 
points at a common, or similar, type of deactivation 
process for both classes of C 2 -symmetric metallocene 
catalysts. This aspect is discussed in part 2 of this work, 
where the formation and reactivation of Zr(allyl) species 
is analyzed. 

Proton NMR analysis of the unsaturated end groups 
of the j-PP samples prepared with 1/MAO in the 
presence of hydrogen shows that the main intramolecu- 
lar chain transfer reaction is /J-CH3 elimination, while 
the microstructural analysis of the polymer samples by 
13 C NMR spectroscopy revealed an unexpected result 
on the composition of the saturated end groups. As 
shown in Figure 1 and Table 2, along with the peaks of 
the "normal" chain-end groups (i.e., the />propyl group, 
A, isobutyl end group, B, and traces of allyl end group, 
D), a new set of signals (C) appeared in the 13 C NMR 
spectra of polypropene samples. 

From the quantitative analysis of 13 C NMR spectra 
(see Table 2) we can observe that, as expected, the 
amount of isobutyl and /?-propyl end groups increases 
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Scheme 1 
p-CH 3 transfer 




Zr— CH 3 + 




D 



terminating 
chain-end 



B 



with hydrogen pressure. The same behavior can be 
noticed in the case of the new group (C in Table 2). 
Furthermore, the A/C ratio is almost constant (^5) in 
samples 3-5. These facts suggest that the new group 
is also a chain-end. 

To better understand the nature of this group, we 
have to make a balance between the initiating and 
terminating chain ends. As the chain-end B can be 
formed in two ways (see Scheme 1), i.e., from ^-methyl 
transfer (it is an initiating chain-end) or from hydrogen 
transfer (it is a terminating chain-end), we have to split 
the two different contributions as [B te rm] — [B] - [BiniJ 
= [B] - [D] (where B te rm refers to the fraction of isobutyl 
group from hydrogen transfer, and Bmit is the fraction 
from reinitiation after /?-methyl transfer). 

As the relation C + A^B-Dis verified for samples 
2-6 (see Table 2), we can safely assign group C to an 
initiating chain-end associated with the direct insertion 
of propene in a Zr-H bond. The complete identification 
of this new chain- initiating group was performed using 
monodimensional (DEPT) and bidimensional (GRASP- 
HSQC and GRASP-HMBC) NMR techniques. 7 The 
more reliable structure for this chain-end group derives 
from a regioirregular propene insertion in the Zr-H 
bond followed by regio- and stereospecific propene 
insertions (see below) according to Scheme 2. 

The same set of signals was observed by Busico and 
Chadwick in 13 C NMR spectra of some polypropenes 
obtained from heterogeneous systems. 8 Independently, 
Randall reported on a similar assignment scheme in the 
case of propene polymerization by using less regio- 
specific metallocenes in the presence of hydrogen. 9 

(7) (a) Bendall, M. R; Doddrell, D. M.; Pegg, D. T. J. Am. Chew. 
Soc. 1981, 103, 4603. (b) Schenker, K. V.: v. Philipsborn. W. J. Magn. 
Reson. 1986. 66, 219. (c) Kay, L. E.; Keifer, P.; Saarinen. T. J. Am. 
Chem. Soc. 1992, 114, 10663. (d) Kay, L. E.; Keifer, P.; Saarinen, T. J. 
Am. Chem. Soc. 1992, 114, 10663. 

(8) Busico, V.; Chadwick, J. C. Personal communication. 
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Table 3. 13 C and *H NMR Assignments 

chemical shift (ppm) 



position* 


type 


carbon 


hydrogen 
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CH 3 


17.76 


0.89 
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CH 


31.93 


1.67 
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CH 3 


20.56 


0.96 


4 


CH 3 


16.30 


0.89 
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CH 


36.47 


1.51 


6 


CH 2 


43.05 


0.93/1.42 



a See atom numbering in Scheme 2 

The complete *H and 13 C NMR assignments are 
reported in Table 3. The chemical shift difference of the 
two diastereotopic hydrogen atoms on carbon 6 (Ad = 
0.49 ppm) can help us to assign the relative configura- 
tion of the two nearest methines. In fact, by comparing 
our value with the values found by Mizuno 10 for a 
similar CH 2 group next to a 2,1-erythro (CH 2 in a meso 
dyad) or a 2,1-threo (CH 2 in a racemic dyad) regio- 
irregular propene unit (see Scheme 3), it is evident that 
in group C the CH 2 is in a meso dyad. This indicates a 
highly stereospecific propene insertion after the initial 
regio mistake. 

With the aim of rationalizing the above experimental 
evidence that a large fraction (about 20%) of propene 

(9) Randall, J. C; Ruff, C. J.; Vizzini, J. C; Speca. A. N.; Burkhardt, 
T. J. In Metalorganic Catalysts for Synthesis and Polymerisation-, 
Kaminsky, W., Ed.; Springer-Verlag: Berlin, 1999; p 601. 

(10) Mizuno. A., Tsutsui. T.; Kashiwa. N. Polymer 1992, 33, 254. 
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insertions on the [rac-Me2C(3-r-Bu-l-Ind)2Zr-H] + spe- 
cies proceeds in a secondary fashion, a DFT quantum 
mechanics/molecular mechanics study (see computa- 
tional details) on the evolution of the analogous [H2C(3- 
r-Bu-l-Ind)2Zr-H] + /propene system has been per- 
formed. 11 This approach has been successfully employed 
in a number of theoretical studies and is based on the 
assumption that the core of the molecular system (that 
is the active site, at which processes of bond breaking 
and forming take place) must be described by an 
electronic structure calculation, while the rest of the 
system (perturbating the active site only through steric 
hindrance) can be described by the simpler (and com- 
putationally faster) molecular mechanics treatment. 

Theoretical Investigation on the Evolution of 
the Zr-H Species. A metal-hydride bond is generated 
upon chain-transfer reaction with molecular hydrogen. 
This highly unsatured species has a strong tendency to 
react with a new monomer molecule: test calculations 
show that the olefin in the [L^ZrH (propene) ] + complex 
(where L 2 = (i?,/^-H 2 C(3-f-Bu-l-Ind)2) has a coordina- 
tion energy of nearly 30 kcal/mol. We are aware that 
our model is largely approximated: in fact, the system 
is considered in vacuo, and moreover, the presence of a 
counterion, as well as the influence of entropic factors, 
is neglected. However, Ziegler and co-workers have 
recently given good reasons for neglecting the influence 
of the counterion and also reported that the TAS value 
for ethylene uptake is about -10 kcal/mol. 12ab There- 
fore, we can reasonably assume that coordination of a 
propene molecule on the [L2Zr-H] + species is a largely 
favored and possibly nonactivated process, even though 
its exothermicity could be overestimated. 

In principle, a propene molecule can coordinate to the 
Zr-H species in four ways: primary (1 ,2) and secondary 
(2,1), both with either re or si enantiofaces. The opti- 
mized structures of the [L^Zr— H (propene)]" 1 " species in 
which the monomer molecule is 1,2-re (model la) and 



(11) Computations have been performed on t/?,/3-H2C(3-f-Bu-Ind)2 
for simplicity. Actually, rac-H2C(3-f-Bu-Ind)2ZrCl2 is a real and 
improved version of 1 that maintains all features typical of 1/MAO. 
Indeed, test calculations with the (/?,i?)-Me2C(3-f-Bu-Ind)2 ligand gave 
quite comparable results; for details, see: Resconi, L.; Balboni, D.; 
Dang, V. A.; Yu, L. C. PCT WO 98/43989. Resconi, L.; Balboni, D.; 
Baruzzi, G.; Fiori, C; Guidotti, S. Submitted to Organometallics. 

(12) Deng, L.; Ziegler, T.; Woo, T. K.; Margl, P.; Fan, L. Organo- 
metallics 1998, 17, 3240. (b) Moscardi G.; Woo, T. K.; Deng, L.; Ziegler, 
T. A Combined DFT QM/MM Study on the Associative Displacement 
of the Terminated Polymer Chain in Ni(II) Diimine-Catalyzed Ethylene 
Polymerization. Unpublished work. 
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1,2-sicoordinatec^^^pel lb) are shown in Figure 2. In 
the same figure, mStm 1c shows the optimized structure 
of [L2Zr-H(propene)] + with a coordinated 2,1-si pro- 
pene. All attempts at optimizing the structure including 
a 2,1 -re propene inexorably led to the 2,1 insertion 
product [L 2 Zr (isopropyl)] 4 " (model 2b in Figure 2), As a 
possible rationalization of this behavior, we recall that 
the enantioselectivity of the regioirregular propene 
insertions into a growing polypropene chain was ex- 
plained in terms of nonbonded interactions between the 
ligand and the methyl group of propene. 13 It was shown 
that in models with an (R,R)-ansa-ji-ligand, the coor- 
dination of a secondary 2, 1 -re monomer is disfavored by 
repulsive interaction with the ligand (see Figure 4 in 
ref 13): similarly, geometries in which the growing 
chain is replaced by a H atom show small nonbonding 
distances between the methyl group of a coordinating 
2,1 -re propene and one six-membered ring of the ligand. 
As a consequence, the olefin is "pushed" to insert into 
the Zr-H bond, giving rise to the more compact 
Zr(isopropyl) insertion product 2b. In Table 4, a splitting 
of the energy differences in terms of their quantum 
mechanics and molecular mechanics components for all 
the species considered in this work is reported. For the 
Zr-H (propene) complexes, the most stable structure 
corresponds to the secondary si propene coordination 
(model lc). This effect is electronic in nature and is due 
to an agostic bond between the metal and one of the 
hydrogen atoms of the methyl group of propene, whereas 
this interaction is not allowed in models (la and lb) 
for primary monomer coordinations. The occurrence of 
this agostic interaction in models including a propene 
molecule 2,1 coordinated to Zr-H was first reported by 
Prosenc and Brintzinger in a modelistic study on the 
isomerization process of the n-propyl growing chain on 
the [Cp 2 ZrC 3 H 7 ] + species. 14 

Insertion into the Zr-H Bond. Some authors have 
already investigated the process of hydrogen transfer 
to the metal center from /7-propyl units 14 " 16 and the 
insertion of coordinated 2,1 propene into the Zr-H bond 
to give an isopropyl unit. 1416 Since there is a substantial 
agreement between the energy barriers shown in the 
mentioned studies (particularly with those of ref 14) and 
our results, we just summarize our findings in Table 4 
and recall that both primary and secondary insertions 
occur with small activation barriers (less than 1 kcal/ 
mol). In the Zr (isopropyl) complex, both the isopropyl 
methyls show an agostic interaction with the zirconium 
and are located at nearly the same distance to the metal 
center (see different views of 2b in Figures 2 and 7). 17 
This conformation provides the best interaction with the 
low-lying metal-centered LUMOs localized in the "equa- 
torial" plane between the Cp ligands. 18 Figure 3 depicts 
the relative energy along the two possible "primary" and 
"secondary" pathways representing the evolution of the 
[L 2 Zr-H]+/propene system: since the reactive Zr-H 
species always reacts with the monomer with small 
activation barriers, primary and secondary propene 

(13) Guerra, G.; Cavallo, L.; Moscardi, G.; Vacatello, M.; Corradini, 
P. J. Am. Chem. Soc. 1994, 116, 2988. 

(14) Prosenc, M. H.; Brintzinger, H. H. Organometallics 1997, 16, 
3889 

(15) Yoshida, T.; Koga. N.; Morokuma. K. Organometallics 1995, 
14, 746. 

(16) Thorsaug, K.; Stovneng J. A.; Rytter, E.; Ystenes, M. Macro- 
molecules 1998, 31, 7149. 
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Figure 2. Models for propene coordination on the [L 2 Zr-H] + complex in which the aromatic ligand is (R,R) coordinated 
to the metal and related insertion products. Hydrogen atoms on the ligand framework have been omitted for clarity. 
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Table 4. Relative Energies (kcal/mol) of the 
Species Involved in the Studied Processes 



A£"(Tot) 



A£(QM) 



A£(MM) 



Zr-H(propene) Complexes 3 
lc 0.0 0.0 

la 3.4 3.4 

lb 3.8 3.6 

Propene Insertion into Zr-H a 
TS[la-2a] +3.5 +3.4 

2a -5.2 -5.4 

TS[lc-2b] 0.3 -0.2 

2b -8.1 -8.2 



-8.1 



Figure 3. Reaction between cationic hydride complex and 
propene leading to the formation of /7-propyl and isopropyl 
units. Solid and dotted lines refer to processes in which si 
and re propene are involved, respectively. Insertion barriers 
are not shown. 

insertions should have the same probability. This result 
does not agree with our experimental data, which 
clearly indicate that the majority (^80%) of the initiat- 
ing groups are /7-propyl (Table 2). In the following, we 

(17) As a main difference with calculations of refs 14 and 16, we 
find that the 2,1 insertion product— the Zr(isopropyl) species 2b— is 
more stable than Zr(/7-propyl) 2a. To check to what extent our findings 
depend on both the particular ligand and the computational technique, 
we compared the relative stability and the geometric features of the 
[Ligand Zr(/?-propyl)] + and [Ligand Zr(isopropyl)] + complexes (Ligand 
= Cp 2 , H 2 C(Cp) 2 , rac-[H 2 C(l-Ind) 2 ]): while for Ligand = Cp 2 we 
reproduced the results shown in refs 14 and 16; with all the bridged 
ligands we again found a higher stability of the Zr(isopropyl) species. 
In our opinion, this is due to the more open ligand geometry induced 
by the bridge: in fact, while for the Cp 2 the angle between the mean 
planes of the five-membered rings is about 50°, for all the bridged 
ligands this value is more than 70°. Evidently, repulsive steric 
interactions in the unbridged catalysts disfavor the more electronically 
stable Zr(isopropyl) species. 

(18) Lauher, J. W.; Hoffmann, R. J. Am. Chem. Soc. 1976, 98, 1729. 



Propene Insertion into Zr(n-propyl) 6 
3a -10.3 -10.9 

max[3a— 3b] -6.1 -7.0 

3b -10.1 -10.5 

4a -9.8 -9.5 

insertion product -20.5 -19.0 

Propene Insertion into Zr (isopropyl) c 
TS[2b-2c] 9.0 8.3 

2c 5.9 5.7 

3c -5.3 -6.6 
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insertion product -15.9 -14.4 
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Propene Rotation in the Two-Step Isomerization 3 

8.1 8.6 -0.5 

8.0 9.2 -1.2 



Displacement of 2,1 Propene from lc a 
-8.3 -7.9 



-0.4 



a Energy of lc chosen as reference value. b Energy of 2a chosen 
as reference value. c Energy of 2b chosen as reference value. 

will show that the way the isopropyl interacts with the 
metal has a detrimental effect on the propagation 
reaction and will provide a rationalization for the 
"dormancy" of the catalytic site bearing an isopropyl 
group. 

Comparative Study on the Propagation Step for 
Catalytic Sites Bearing a u-Propyl and Isopropyl 
Chain. Propene Insertion into Zr(n-propyl). The 
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Figure 4. Concerted propene coordination- chain rear- 
rangement process for the formation of the preinsertion 
intermediate (3b) starting from the olefin adduct (3a) and 
transition state for propene insertion into the Zr— C3H7 
bond (4a). The hydrogen atoms on the ligand framework 
have been omitted for clarity. Thicker lines represent bonds 
above the plane of the paper. 
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Figure 5. Energy profile for the concerted propene 
coordination-chain rearrangement process from 3a to 3b. 

ethylene insertion step onto a linear growing chain has 
been investigated in several theoretical works. 151619 - 20 
In particular, the insertion of ethylene in the presence 
of the bisCp and bis(MesCp) ligands has been among 
the subjects of ref 16, in which the authors also 
investigated the process of chain propagation starting 
from a y-agostic n-propyl. However, other authors 
evidenced that the direct precursor to olefin insertion 
is the chain with a /?-agostic bond, since y-agostic bonds 
are replaced by /S-agostic ones prior to insertion. 20 
Moreover, calculations aimed at optimizing a structure 
in which a si-propene, a y-agostic isobutyl chain, and a 



(19) Lohrenz, J. C. W.; Woo, T. K.; Ziegler, T. J. Am. Chem. Soc. 
1995 117 12793. 

(20) Woo. T. K.; Margl, P. M.; Lohrenz, C W.; Bldchl, P. E.; Ziegler, 
T. J. Am. Chem. Soc. 1996, 118, 13021. 
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Figure 6. Geometric aspects of the concerted chain 
rearrangement— monomer coordination processes leading 
to the formation of the preinsertion intermediate for the 
primary (1,2) monomer insertion. 

(i?,i3-Me 2 C(3-r-Bu-Ind) 2 ligand are bonded to a zirco- 
nium center resulted in the formation of the /?-agostic 
isobutyl, si-propene 7r-complex. 21 Therefore, we chose not 
to further investigate the alternative insertion route 
proceeding from the y-agostic conformation of the n- 
propyl. 

Experimental data show that monomer insertion into 
both /7-propyl and isopropyl initiating groups is primary 
and occurs with the same preferred enantioface of the 
successive insertion steps. Moreover, as already pointed 
out in a recent molecular mechanics study by some of 
us, 22 the si enantioface has been found to be the "right" 
one for (#ig-Me 2 Si(3^Bu-l-Ind) 2 ZrCl 2 . Therefore, we 
have restricted our study to the primary si-propene 
propagation starting from the /J-agostic structure 2a in 
the presence of the sterically crowded ligand framework 
provided by the (i?,^-[H 2 C(3-£-Bu-Ind) 2 ]ZrCl 2 catalyst. 

In the presence of the sterically unhindered bis- 
(cyclopentadienyl) ligand, the formation of a ^-bonded 
ethylene, /3-agostic /7-propyl complex as a precursor for 
chain propagation has been found by Rytter et al. to 
proceed without energy barrier. 16 The same authors 
found that to form an analogous complex in models 
including the crowded bis(pentamethylcyclopentadienyl) 
ligand, an ethylene molecule has to be "pushed* into the 
ligand framework, overcoming an energy barrier of 
nearly 9 kcal/mol. Our attempts to optimize the struc- 
ture of a zirconocene model complex with the [H 2 C(3- 
f-Bu-l-Ihd) 2 ] ligand, a 77-propyl ^-agostic chain, and a 
jr-coordinated si-propene without any geometry con- 
straint led to the formation of an adduct (labeled 3a in 
Figure 4) which is stabilized by nearly 10 kcal/mol with 
respect to the /?-agostic complex 2a plus free propene 
system. This complexation enthalpy is considered suf- 
ficiently exothermic to guarantee a good catalytic activ- 
ity. 23 Then, the rearrangement of the growing chain 
toward the a-agostic conformation and olefin coordina- 
tion to form a preinsertion intermediate (3b in Figure 
4) proceeds in our opinion in a concerted fashion. In 
Figures 5 and 6, the relative energy and the variation 
of several geometrical parameters as Zr-C 2 , Zr-H a , and 

(21) Resconi, L. Fait, A.; Piemontesi. F.; Camurati, I.; Moscardi, G.; 
Cavallo, L.; Mllano, G. Part 2 of this work, manuscript in preparation. 

(22) Toto, M.; Cavallo, L.; Corradini, P.; Moscardi. G.; Resconi, L.; 
Guerra, G. Macromolecules 1998, 31, 3431. 

(23) Margl, P.; Deng, L.; Ziegler, T .Organometallics 1998. 17, 933. 
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propene-bound intermediate 
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Zr-rfy distances (atom labeling: see 3a in Figure 4) and 
the dihedral angle Ci-Zr-Ca^C^ (labeled #i) are de- 
picted versus the Zr-C a -C^ angle (in the following 
labeled as 4>), which has been found to be the best 
reaction coordinate. 16 The /?-agostic bond is almost 
completely lost at O = 100°, but due to repulsive 
interactions between the chain and the monomer, the 
jr-complex cannot be completely established. Moreover, 
the chain conformation is not yet suitable for an 
a-agostic interaction. Only at 3> = 120° is the formation 
of a jr-olefin, a-agostic complex allowed, as one can see 
from the sudden variation of #i and the other relevant 
distances. Further opening of only reduces the steric 
repulsions between monomer and chain, but geometric 
parameters do not change significantly. Starting from 
the preinsertion intermediate 3b, the decrease of the 
C2-C a distance gives rise to the y-agostic insertion 
product. Model 4a in Figure 4 depicts the transition 
state for monomer insertion. Its energy is only 0,30 kcal/ 
mol above that of the preinsertion intermediate 3b. 
Summing up, the /?-agostic — a-agostic rearrangement 
process in concert with monomer coordination is the 
rate-determining step for propagation. The calculated 
energy barrier (4.2 kcal/mol, obtained as the difference 
between the energy level of max[3a-*3b] and that of 
3a in Table 4) is similar to that found in analogous 
calculations performed on the (7?,^)-[Me2C(3-t-Bu-l- 
Ind)2Zr(si-propene) (isobutyl)]" 1 " model. 21 

Propene Insertion into Zr(isopropyl). Propene 
insertion can start from an alkene-bound state in which 
the isopropyl can be either a-agostic or /J-agostic (see 
Scheme 4). In the former, the methyl group of the 
propene is syn (i.e., close) to one methyl of the isopropyl, 
while in the /S-agostic geometry these groups are in an 
anti disposition. A comparison by molecular mechanics 
calculations 22,24 between approximate transition states 
for propene insertion into a primary growing chain in 
which the propene methyl group and the chain are syn 
and anti, respectively, indicated that the syn transition 
state is 3-5 kcal/mol higher in energy. Therefore, we 
assume that the propene insertion into the Zr(isopropyl) 
species with an a-agostic interaction is disfavored and 

(24) Kawamura-Kuribayashi, H.: Koga, N.; Morokuma. K. J. Aw. 
Chem. Soc. 1992, 1 14, 8687. 
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Q = 107.6°* 
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Figure 7. Rearrangement of the isopropyl group from the 
double to the single /?-agostic conformer. Thicker lines 
represent bonds above the plane of the paper. Ligand 
hydrogens are omitted for clarity. Q is the improper torsion 
angle Zr-C^i-C a -C^ 2 . Xcp and Xcp' indicate the centroids 
of the five-membered rings. 

restrict our investigation to the insertion step into the 
/?-agostic conformer of the [LaZrftsopropyl)] 4 * species. 

As remarked before, our models do not represent the 
real system: in particular, the presence of a high 
amount of monomer (polymerization tests were per- 
formed in liquid monomer) is not accounted for. How- 
ever, our goal is not to evaluate the energy stabilization 
of metal complexes due to interactions with surrounding 
monomer (which, in principle, can take place for all the 
considered models), but to verify whether the monomer 
has a role into a specific reaction of the catalytic cycle. 
In the previous section, we have shown that the re- 
arrangement of chain conformation (from fi- to a-agostic) 
and monomer coordination could be concerted processes 
starting from an olefin adduct. The insertion reaction 
into the Zr(isopropyl) bond has been investigated by 
looking for an analogous concerted channel; however, 
the conformation of the isopropyl (as discussed above) 
in the alkene-free state causes the monomer coordina- 
tion step to be quite different. The high stability of the 
[L2Zr(isopropyl)] + species and the strong nonbonded 
interactions between the olefin and the closest methyl 
group of the isopropyl prevent the formation of a 
catalyst-propene adduct; in fact, all attempts to opti- 
mize such a structure without geometry costraints 
always resulted in the rejection of the monomer (Zr- 
Coiefln > 5 A). This suggests that the propagation 
reaction probably proceeds through a rearrangement of 
the isopropyl in order to make room for the incoming 
monomer. Indeed, we investigated such a process by 
choosing as a reaction coordinate the improper torsional 
angle Zr-C^i-Ca-C^ 16 (labeled Q, see model 2b in 
Figure 7 for atom labeling). Some optimized structures 
calculated for the rearrangement step are shown Figure 
7; it is worth noting that initially the isopropyl tends to 
keep its symmetrical conformation. Only at Q ^ 110° 
does a sudden shift toward a different conformation 
occur, as one can see also on inspection of models in 
Figure 7 and from the plots of Figures 8 and 9, in which 
the trends of the relative energy and of some relevant 
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Figure 8. Energy profile for the rearrangement of the Zr- 
(isopropyl) from 2b to 2c (see Figure 7 for atom labeling). 
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Figure 9. Geometric aspects of the unimolecular rear- 
rangement of Zr(isopropyl) from 2b to 2c (see Figure 7 for 
atom labeling). 

geometric parameters are reported for the conversion 
of 2b into 2c. Test calculations aimed at finding a stable 
olefin ^-complex at Q = 100° and 110° confirm that the 
monomer does not take part in the isopropyl rearrange- 
ment process until a conformer with a single /?-agostic 
interaction of the isopropyl (2c, in which # is nearly 
130°) is formed. In fact, we find that at Q = 110° the 
olefin carbons are still at nearly 5 A from the metal. 
The energy profile of Figure 8 shows that the loss of 
the agostic interaction between one of the /^-hydrogens 
with the metal center passes over an energy barrier of 
9.0 kcal/mol. The single /?-agostic species 2c is able to 
coordinate a monomer, leading to the formation of the 
preinsertion intermediate 3c shown in Figure 10. The 
insertion step then follows the known C2-C a reaction 
coordinate, with an activation energy of 4.4 kcal/mol (see 
Table 4 and model 4b in Figure 10), leading to a 
5-agostic insertion product. The propagation into the Zr- 
(isopropyl) bond looks similar to the so-called "backside 
insertion" already discussed in various papers for the 
polymerization of ethylene. 1619 20 In particular, Ziegler 
and co-workers evidenced in ref 19 the influence of 
entropic effects, which can further disfavor this process 
in zirconocene catalysts. Therefore, the activation bar- 
rier for monomer insertion into the Zr(isopropyl) bond 
could be underestimated. 

Ratio between "Initiating" Groups. Figure 11 
summarizes our results on the evolution of the [I^Zr- 
H] + /propene system, proceeding from the first monomer 
insertion into the Zr-H bond to the successive one. By 
comparing the path after the initial 2,1 propene inser- 
tion into the Zr-H bond with the one that follows the 




<SS> 




Figure 10. Preinsertion intermediate 3c and the transi- 
tion state 4b for propene insertion into the Zr(isopropyl) 
bond. Thicker lines represent bonds above the plane of the 
paper. Ligand hydrogens omitted for clarity. 



[LZr-H] + propene 



TS|2b-2c) jTTL 




Figure 11. Schematic reaction profiles for propene inser- 
tion into Zr-H and Zr-alkyl bonds. Energy differences in 
kcal/mol. 

initial 1,2, it is evident that catalytic sites bearing 
/7-propyl and isopropyl groups do have much different 
propagation barriers. 

In particular, the formation of the olefin adduct 3a 
starting from Zr-H and free propene does not have a 
significant activation barrier, further indicating that 1,2 
propene insertion into Zr-H is irreversible. On the 
contrary, the Zr(isopropyl) species 2b must rearrange 
to a less stable conformer 2c before reacting with the 
monomer. The "dormant" nature of 2b can prevent the 
formation of a nearly equal number of "initiating" 
/7-propyl and isopropyl groups (as the results on the 
insertion into Zr-H could suggest) if the "slow" site 
undergoes side reactions. As shown in Scheme 5, in the 
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presence of hydrogen the Zr(isopropyl) species could be 
decomposed in propane and a Zr-H species. Other 
options arise form the evidence that the "slow" site 2b 
can propagate through rearrangement into 2c or un- 
dergo hydride transfer back to lc with comparable 
energy barriers (that is, the 2,1 propene insertion into 
the Zr-H bond is reversible). We now analyze the 
possible mechanisms for the 2b — * 2a isomerization. In 
principle, this can occur in three different ways (exclud- 
ing unimolecular fi-H transfer to Zr and dissociation of 
the 2,1 propene): (1) bimolecular ji-H transfer to the 
monomer; this mechanism can be ruled out by the fact, 
outlined in the previous section, that the rearrangment 
of 2b generates the "backside" monomer coordination 
(that is, propene anti to the /?-agostic bond, see Figures 
7 and 10), which is not a suitable precursor of a 
bimolecular hydrogen exchange process; (2) fi-H transfer 
to Zr, rotation from secondary to primary of the coor- 
dinated propene, then reinsertion into the Zr-H bond; 
(3) /?-H transfer to Zr and associative displacement by 
a new monomer molecule. We have tested these last two 
alternatives. 

Direct Isomerization. The simplest hypothesis for 
a mechanism that changes a 2,1 -coordinated a-olefin 
into a primary (1,2) one is in-plane rotation. 1416 This 
process is part of the two-step mechanism proposed by 
Chien 25 for the 2,1 — 3,1 isomerization of the secondary 
growing chain-end, which has been observed in poor 
regiospecific zirconocenes. 25-29 In this section, we briefly 
investigate the influence of the [H2C(3-f-Bu-l-Ind) 2 ] 
ligand framework, which is both bridged and Cz sym- 
metric, on the energetics of propene rotation that shifts 
the [L2Zr-H(2,l -propene)] + species (lc) into the [L^Zr- 
H(l,2-propene)] + one (lb). In Figure 12, the total QM/ 
MM energy is plotted as a function of the H-Zr-X 0 i- 
Ci dihedral angle (labeled #o; Xoi indicates the centroid 
of the propene double bond, see Figure 13 for atom 
labeling), while in Figure 13 are shown the structures 
5a and 5b, corresponding to the two maximums of the 
rotation path. The strong nonbonded interactions be- 
tween the rotating propene and the ligand force the 
olefin to interact with the metal in a highly unsym- 
metrical fashion. It is our feeling that in this situation 
of high energy the rotating olefin could be easily 

(25) Rieger, B.; Chien, J. Polym. Bull 1989, 21, 159. 

(26) Resconi, L.; Fait, A.; Piemontesi, F.; Colonnesi, M.; Rychlicki, 
K; Zeigler, R. Macromolecules 1995, 28, 6667. 

(27) Busico, V.; Caporaso, L.; Cipullo. R.; Landriani, L.; Angelini, 
G.; Margonelli, A.; Segre, A. L. J. Am. Chem. Soc. 1996, 118, 2105. 

(28) Soga, K.; Shiono, T.; Takemura, S.; Kaminsky, W. Makromol 
Chem., Rapid Commun. 1987, 8, 305. 

(29) Schupfner, G.; Kaminsky, W. J. Mol Catal A: Chem. 1995, 
102, 59. 
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Figure 12. Energy profile for the step of propene rotation 
(with respect to the H-Zr-Xoi-C2 dihedral angle) in the 
two-step isomerization mechanism. 




Figure 13. Selected geometries corresponding to points 
5a and 5b of Figure 12. Thicker lines represent bonds 
above the plane of the paper. Ligand hydrogens are omitted 
for clarity. 

displaced by a new monomer molecule. To our knowl- 
edge, only Prosenc and Brintzinger have so far proposed 
a mechanism that is alternative to in-plane olefin 
rotation. 14 

Associative Displacement. The reaction between 
the [L 2 Zr-H(2,l-propene)] + species (lc) and a primary 
propene could either reverse the hydrogen-transfer 
process back to the "dormant" species [L2Zr(isopropyl)] + 
(2b) or, alternatively, generate the [L 2 Zr(/>propyl)] + 
complex upon primary propene insertion into the Zr— H 
bond with concurrent dissociation of the secondary 
propene, that is, insertion by associative displacement 
of the 2,1 propene. 30 This latter hypothesis is supported 

(30) Johnson, L. K.; Killian, C. M.; Brookhart, M. J. Am. Chem. Soc. 
1995, 117, 6414. 
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•r olefin dissociation is about 10 



Figure 14. Reaction product between free propene and 
complex lc in the displacement of the 2,1 -coordinated 
monomer. Ligand hydrogens are omitted for clarity. 

by the extremely high regiospecificity of the catalyst, 
which suggests that in the presence of a growing chain 
the coordination of a 2,1 propene is highly disfavored. 
In other words, by this sequence of events, a bimolecular 
self-correcting isopropyl to /7-propyl isomerization pro- 
cess occurs. To determine the energetic cost of the 
approach of the 1,2 propene to the inner coordination 
sphere of the metal, we tried to optimize a structure 
for a double olefin metal-hydride intermediate in which 
the Zr-H bond is placed between a 1,2-si and a 2,1-si 
propene molecule. Every attempt always resulted in the 
formation of the [L2Zr(/7-propyl)(2,l-propene)] + complex 
(model 6 in Figure 14), Moreover, the exothermicity 
(13.5 kcal/mol) of the process from lc plus free propene 
to 6 is higher than the TAS at room temperature for 
olefin uptake (about 10 kcal/mol, l2a b ). The above results 
indicate that the overall process lc + 1,2-propene — * 6 
could have a fairly small activation energy. In complex 
6, the 2,1 propene is loosely bonded, and most impor- 
tantly, there is no chance for its insertion into the Zr- 
(/7-propyl) bond. We evaluate that in model 6 the 
coordination energy of the 2,1 propene is nearly 8 kcal/ 
mol: therefore, the complete dissociation of the 2,1 
propene may be a spontaneous process, since the TAS 



at room temperatj 
kcal/mol. 12ab 

The results of the two proposed schemes for the 
conversion of 2,1 units into 1,2 ones are summarized in 
Figure 15: the energy barrier for the two-step isomer- 
ization (nearly 16 kcal/mol) starting from the double 
/9-agostic species 2b indicates that this process could be 
largely disfavored and that conversion of "slow" into 
"fast* sites (Scheme 5) could be due to hydrogen or, 
alternatively, to a process that reverses the initial 2,1 
propene insertion, followed by the displacement of the 
2,1 propene by a new 1,2 coordinating monomer. 

That a self-correcting isopropyl to /7-propyl isomer- 
ization process can also occur without hydrogenolysis 
is indicated also by the fact that the 2,3-dimethylbutyl 
group is never observed in i-PP samples prepared with 
1/MAO without hydrogen, even at the lowest propene 
concentrations, where Zr-H bonds are generated by 
unimolecular chain transfer to the metal. 6b ~ d As a final 
remark, the substantial equivalence between the energy 
barrier for chain propagation and for 2,1 propene 
elimination from lc (Figure 15) suggests a [1,2]/[2,1] 
ratio of about 3, in the absence of hydrogenolysis. 

Conclusions 

Addition of molecular hydrogen to the polymerization 
of liquid propene catalyzed by rac- [Me2C (3- f-Bu-l-Ind) - 
ZrC^j/MAO produces a relatively small increase in 
catalytic activity, while the polymer molecular weight 
is strongly reduced: polypropene number average mo- 
lecular weight decreases from about 40 000 without 
hydrogen to about 2000 when P$\i) = 4 bar. 

A new saturated chain-end group is observed in the 
13 C NMR spectra of these samples, which we assign to 
secondary propene insertion into the Zr— H bond, gener- 
ated after chain transfer to hydrogen. The ratio between 
primary to secondary initial insertions ([1,2]/[2,1]) was 
found to be nearly equal to 5 and constant over the 
investigated range of hydrogen partial pressure. 

After the initial insertion the catalyst becomes re- 
giospecific again, as no 2,1 units inside the polymer 
chain are detectable in the 13 C NMR spectra. 

Theoretical calculations provide an explanation for 
the observed ratio between "initiating" polypropene end 
groups in terms of the lower reactivity of the Zr- 
(isopropyl) site compared to that of the Zr(/7-propyl) 
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Figure 15. Schematic representation of the proposed mechanisms for the conversion of the "slow" Zr(isopropyl) site into 
the "fast" Zr(/>propyl) one. Relative energies in kcal/mol. 
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one: the less reactive Zr(isopi|^Hp undergoes competing 
chain-transfer reactions. As ^^theoretical results on 
the evolution of the [L 2 Zr-H] + in the presence of 
propene have been obtained for a highly crowded model 
catalytic system, we believe that the same reactions can 
occur in less hindered catalysts, and therefore this 
behavior should not be limited to the present catalytic 
system. Indeed, experimental observations are accumu- 
lating that show the occurrence of 2, 1 propene insertion 
into the Zr-H bond for both homogeneous 9 and hetero- 
geneous 8 Ziegler-Natta catalysts. 

The "dormancy" of the [L 2 Zr(isopropyl)] + can be 
rationalized by the high stability arising from a double 
/?-agostic interaction. It is now straightforward to sug- 
gest that a similar effect can occur in the propagation 
into a secondary growing polypropene chain (which can 
be modeled with a sec-butyl group), and this would also 
explain the decrease in the propagation rate and the 
occurrence of chain-transfer reactions observed after an 
occasional regioirregular insertion in poor regiospecifk 
catalysts. 

The most likely Zr(isopropyl) to Zr(n-propyl) inter- 
conversion pathway includes (1) unimolecular fi-H 
transfer to the Zr atom with formation of a Zr-H(2,l- 
propene) cation, which is stabilized by an agostic 
interaction between one of the methyl protons and the 
metal center, and (2) associative displacement with an 
incoming primary propene molecule. 

As a final remark, the existence of a concerted chain 
conformation rearrangement-monomer coordination 
channel leading to primary propene insertion into a 
primary growing chain can contribute to rationalizing 
the influence of ligand structure on processes of chain 
transfer and therefore on the molecular weight of the 
polymer; in fact, the formation of a ^-propene, /J-agostic 
chain complex as a precursor of the hydrogen exchange 
reaction seems to be an unfavored process for highly 
crowded catalysts. As a consequence, other chain- 
transfer reactions, possibly with a unimolecular course, 
may occur. Indeed, H2C(3-f-Bu-l-Ind)2ZrCl2, Me2C(3-t- 
Bu-l-Ind)2ZrCl 2 , and (Me5Cp) 2 ZrCl 2 undergo chain trans- 
fer mainly via /?-methyl elimination, 6 31 whereas less 
crowded catalysts as (Cp) 2 ZrCl2, rac-[Me 2 C(l-Ind) 2 ]- 
ZrCl 2 , rac-[Me 2 C(3-Me-l-Ind) 2 ]ZrCl 2 , and rac-[Me 2 C(3- 
Me3Si-l-Ind)2]ZrCl 2 in combination with MAO termi- 
nate mainly via /?-H elimination. 6 - 31 Moreover, by 
increasing the size of the substituent in the position 3 
of indene, the molecular weights of polymers produced 
with fully regiospecific catalysts of the type rac-[MezC- 
(3-R-l-Ind) 2 ]ZrCl 2 /MAO (R = Me, Me 3 Si-, f-Bu) are 
increased. These experimental data could be rational- 
ized by the lower availability of the monomer to chain 
transfer in highly hindered catalysts. 

Experimental Part 

The catalyst precursor rac-Me 2 C(3-rert-butyl-l-indenyl)2ZrCl 2 
was prepared by Dr. I. Nifantev according to the procedure 
described in ref 6a. 

MAO was a commercial 10% solution in toluene (Witco). The 
zirconocene was prereacted with MAO in toluene for 10 min 
and pressurized into the polymerization reactor containing 
propene. The polymerization tests were carried out as fol- 



(31) Resconi, L.; Piemontesi, F.; Franciscono, G.; Abis, L.; Fiorani, 
T. J. Am. Chem. Soc. 1992, 114, 1025. 
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lows: 0.4 L of pron^^^ was charged into the 1 L stainless 
steel autoclave (Bucl^^30 °C, the temperature of the reactor 
was brought to 50 °C, and, in quick sequence, the catalyst/ 
cocatalyst mixture was loaded under a nitrogen overpressure, 
excess pressure was released, and the desired amount of 
hydrogen (measured into a stainless steel 100 mL vial) was 
added; the polymerization test was carried out for 1 h at 50 
°C and then stopped with CO. The polymers were recovered 
by evaporating the solvents under vacuum and extracted with 
refluxing toluene under nitrogen in a Kumagawa extractor to 
remove catalyst residues, toluene was evaporated, and the 
polymer was dried at 60 °C under high vacuum. 

Proton and carbon spectra were obtained using a Varian 
300 spectrometer operating in the Fourier transform mode at 
120 °C. The samples were dissolved in C2D2CI4. As reference 
the residual peak of C2DHCI4 in the l H spectra (5.95 ppm) and 
the peak of the mmmm pentad in the 13 C spectra (21.8 ppm) 
were used. The carbon spectra were acquired with a 90° pulse 
and 12 s of delay between pulses. About 3000 transients were 
stored for each spectrum. 

Models and Computational Details. Elements of Chiral- 
ity. The elements of chirality, which are relevant for the 
present study, are briefly recalled in order to indicate the used 
terminology. First of all, upon coordination, the prochiral 
propene gives rise to nonsuperimposable si and re coordina- 
tions. 32 Another element of chirality is that of the catalytic 
site, arising from coordinated ligands other than the alkene, 
the monomer, and the growing chain. For metallocenes with 
prochiral ligands we used the notation (hQ or (5), in paren- 
theses, according to the Cahn-Ingold-Prelog rules 33,34 ex- 
tended by Schlogl. 35 Without loss of generality, all the calcu- 
lations reported here refer to models with the {R,R) chirality 
of coordination of the H 2 C(3-tert-butyl-l-indenyl) 2 ligand. 

Computational Details. Stationary points on the potential 
energy surface were calculated with the Amsterdam Density 
Functional (ADF) program system, 36 developed by Baerends 
et al. 37,38 The electronic configuration of the molecular systems 
were described by a triple-^ basis set on zirconium for 4s, 4p, 
4d, 5s, and 5p. Double-^ STO basis sets were used for carbon 
(2s,2p) and hydrogen (Is), augmented with a single 3d and 2p 
function, respectively. 3940 The inner shells on zirconium 
(including 3d) and carbon (Is) were treated within the frozen 
core approximation. A set of auxiliary 41 s, p, d, f, and g STO 
functions, centered on all nuclei, was used in order to fit the 
molecular density and present Coulomb and exchange poten- 
tials accurately in each SCF cycle. Energetics and geometries 
were evaluated by using the local exchange -correlation 
potential by Vosko et al., 42 augmented in a self-consistent 
manner with Becke's 43 exchange gradient correction and 
Perdew's 44 45 correlation gradient correction. First-order scalar 
relativistic corrections were added to total energy. 46 - 47 A recent 
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2,1 Propene Insertion in the Zr-H Bond 
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study by Jensen and Borve on Ti ethylene insertion reactions 
concluded that the functional used in the present paper is well 
suited to reproduce high-level ab initio coupled cluster calcula- 
tions for this type of reactions. 48 

The ADF program was modified 49 " 51 to include standard 
molecular mechanics force fields in such a way that the QM 
and MM parts are coupled self-consistently, according to the 
method prescribed by Morokuma and Maseras. 52 The model 
QM system and the real QM/MM systems are reported in 
Chart 1 . The partitioning of the systems into QM and MM 
parts only involves the skeleton of the catalyst's ligand; that 
is, either the n-propyl, the isopropyl, or the monomer is always 
totally composed by pure QM atoms. Moreover, the metal 
atom, the hydrogen bonded to it (in models for \LiZr- 
H (monomer) ] + ), the five-membered Cp rings of the catalysts, 
and the CH2 bridge connecting the two Cp rings are also 
composed by pure QM atoms. The only MM atoms, hence, are 
the carbon and hydrogen atoms needed to transform the pure 
QM H 2 C(Cp) 2 skeleton of the ligand into the CH 2 (3-te/*-butyl- 
l-indenyl)2 ligand. The connection between the QM and MM 
parts occurs by means of the so-called capping "dummy" 
hydrogen atoms, which are present in the model system only. 52 
These capping atoms are replaced in the real system by the 
corresponding "linking" carbon atom. 52 

The QM and MM parts are thus linked by the "capping" 
hydrogen atoms and coupled by van der Waals interactions. 
The geometry optimization of the whole system was carried 



(47) Snijders, J. G.; Baerends, E. J.; Ros, P. Mol. Phys. 1979, 38, 
1909. 

(48) Jensen, V. R.; Borve, K. J. J. Comput. Chem. 1998, 19, 947. 

(49) Deng, L.; Woo, T. K.; Cavallo, L.; Margl, P. M.; Ziegler, T. J. 
Aw. Chem. Soc. 1997, 119, 6177. 

(50) Woo, T. K.; Cavallo, L.; Ziegler, T. Theor. Chew. Acc. 1998, 100, 
307. 

(51) Cavallo, L.; Woo, T. K.; Ziegler, T. Can. J. Chew. 1998, 76, 1457. 

(52) Maseras, F.; Morokuma, K. J. Cowput. Chew. 1995, 16, 1170. 



Organometallics, Vol 18, No. 25, 1999 5275 

out within this couc^^kcheme between QM and MM atoms. 
In the optimization^We MM part, the Ccp-Cbenzene (respec- 
tively, Ccp-C f -Bu) bonds crossing the QM-MM border were 
kept 0.33 (respectively, 0.41) A longer than the optimized 
Ccp-H distance. This corresponds to optimized C-C^sp 2 ) and 
C-C(sp 3 ) bond distances roughly equal to 1.43 and 1.51 A, 
respectively. Further details on the methodology can be found 
in previous papers. 49,53 

The molecular mechanics potential was described by the 
AMBER95 force field; 54 the carbons of the indene system were 
assigned with atom type "CA", the aromatic ring hydrogen 
atoms with "HA", and the remaining carbon and hydrogen 
atoms with "CT" and "HC", respectively. For the atoms that 
are involved in reactive events at the active site, we followed 
the convention that atoms whose type changes during a 
reaction preserve the original type also in the transition state. 
For example, the propene sp 2 C atoms have a C sp 2 "CA" van 
der Waals parameter both in the reactant and in the transition 
state and an sp 3 "CT" one in the products. Electrostatic 
interactions were not included in the molecular mechanics 
potential. 

All the structures are stationary points on the combined 
QM/MM potential surface, and no symmetry constraints were 
used. Geometry optimizations were terminated if the largest 
component of the Cartesian gradient was smaller than 0.002 
au. AH reported linear transit calculations involve full geom- 
etry optimization along a reaction coordinate that is con- 
strained in each step. Transition state geometries (labeled TS) 
were approached from the reactant side by a linear-transit 
procedure, performing full transition state searches starting 
from geometries corresponding to maxima along the linear- 
transit curves. At the end of each transition state search, the 
approximated Hessian presented one negative eigenvalue. 
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